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We have considered the nonlinear response of mesoscopic systems of non-interacting electrons to 
the time-dependent external field. In this consideration the inelastic processes have been neglected 
and the electron thermalization occurs due to the electron exchange with the reservoirs. We have 
demonstrated that the diagrammatic technique based on the method of analytical continuation or 
on the Keldysh formalism is capable to describe the heating automatically. The corresponding 
diagrams contain a novel element, the loose diffuson. We have shown the equivalence of such a 
diagrammatic technique to the solution to the kinetic equation for the electron energy distribution 
function. We have identified two classes of problems with different behavior under ac pumping. In 
one class of problems (persistent current fluctuations, Kubo conductance) the observable depends 
on the electron energy distribution renormalized by heating. In another class of problems (Landauer 
conductance) the observable is insensitive to heating and depends on the temperature of electron 
reservoirs. As examples of such problems we have considered in detail the persistent current fluctu- 
ations under ac pumping and two types of conductance measurements (Landauer conductance and 
Kubo conductance) that behave differently under ac pumping. 
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I. INTRODUCTION 



Recently there has been a considerable interest in non- 
equilibrium mesoscopics. The effect of adiabatic charge 
pumping [EJ has been experimentally observed H and 
discussed theoretically ||. Weak localization in a quan- 
tum dot under ac pumping has been theoretically stud- 
ied Q. The non-equilibrium noise has been suggested 
[^],[| as a cause of both the low temperature dephasing 
saturation and the anomalously large ensemble aver- 
aged persistent curr ent g. The results of Ref. [| arc 
based on the earlier works M on the ensemble averaged 
dc current caused by the quantum Aharonov-Bohm rec- 
tification of the external ac electric field. Without the 
Aharonov-Bohm magnetic flux the rectified dc current or 
voltage ('photovoltaic effect') has zero ensemble average 
but can exist in individual mesoscopic samples because 
of the specific arrangement of impurities or irregularities 
in the dot's shape. This effect was suggested long ago 
fiof and reconsidered very recently for the case of the 
quantum dot 

Theoretical description of all the effects listed above re- 
quires to go beyond the linear response theory and to con- 
sider the essentially non-linear response to the ac pump 
field. This raises a question on the 'minimal model' for 
the adequate description of nonlinear responses in meso- 
scopic systems. For the linear conductance the minimal 
model is the system of non-interacting electrons with 
the impurity scattering and interaction with the exter- 



nal electric field. Such a model does not explicitly con- 
tain dissipation. Yet it allows to obtain a correct value 
of conductivity that is the key quantity for the dissipa- 
tion function. We will refer to a description based on a 
model of that sort as the 'dynamical approach'. In this 
approach the electron-electron (e-e) and electron-phonon 
(e-ph) interaction is neglected and the stationary regime 
under external pumping is reached in an open system 
via the escape of 'hot' electrons into the massive leads 
playing a role of an electron bath. 

The question we address in this paper is this: To what 
extent this model applies to the nonlinear phenomena in 
mesoscopic systems and how one can see its limitations 
through intrinsic inconsistencies and physical paradoxes. 

Another important issue we address in this paper is 
how to describe heating effects by the impurity diagam- 
matic technique without explicitly solving the kinetic 
equation. It turns out that heating can be described au- 
tomatically by the new class of diagrams containing the 
'loose diffusons' with one free end. They are contrasted 
to the ordinary diffusons and cooperons which are con- 
nected in loops by the 'Hikami boxes' and describe the 
effect of electron phase-coherence. Thus we show the way 
to separate the heating and the dephasing effects on the 
level of the impurity diagrammatic technique. 

The paper is organized as follows. In Sec. II and Sec. Ill 
we discuss the general structure of the perturbation the- 
ory in the external ac field using the method of analyti- 
cal continuation and the Keldysh technique and describe 
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simple rules of the impurity diagrammatic technique in 
the time domain. In Sec. Ill we derive the diffusion prop- 
agators ('diffusons' and 'cooperons') in the external ac 
field at different boundary conditions. Sec. IV is central 
for the paper. There we introduce the 'loose diffusons' 
and demonstrate that evaluating the diagrams with the 
loose diffusons is equivalent to the solution of the kinetic 
equation for the electron energy distribution. We also 
discuss the paradoxes connected with the loose diffusons 
in closed electronic systems. As an example of the role 
of the loose diffusons we consider in Sec.V the variance 
of the persistent current fluctuations under ac pumping. 
For mesoscopic rings connected to an electron reservoir 
by a passive lead we compute the temperature depen- 
dence of the persistent current fluctuations in equilibrium 
and under the harmonic ac pumping. In Sec. VI we con- 
sider the problem of dc conductance under ac pumping 
in two different experimental geometries that correspond 
to measurements of the Landauer and the Kubo conduc- 
tances. We re-derive the expression for the Landauer 
conductance in terms of the electron Green's function in 
the time domain and show that the loose diffusons can- 
not be build in this problem. It means that in systems 
of non-interacting electrons the Landauer conductance is 
insensitive to the electron energy distribution inside the 
mesoscopic system and thus is insensitive to heating. In 
contrast to that the Kubo conductance is sensitive to the 
non-equilibrium electron distribution in the correspond- 
ing system. In Conclusion we summarize the main results 
of the paper and point out on its obvious extensions. 



II. ANALYTICAL STRUCTURE OF THE 
NONLINEAR DYNAMICAL RESPONSE 



In this section we describe the general analytical struc- 
ture of the non-linear response of an arbitrary order 
in the external field using the formalism of the ana- 
lytical continuation jl2j and the Keldysh diagrammatic 
technique |l^,|I^|. We will show that causality encoded 
in the triangular matrix structure of the Green's func- 
tions in the Keldysh technique, allows at most one point 
where the string of retarded electron Green's functions is 
switched to the string of advanced Green's functions in 
the expression for the non-linear response of an arbitrary 
order. 



A. Causality of a non-linear response and the 
method of analytical continuation 



The formalism of analytical continuation is based 
on the explicit assumption of causality. One starts 
with the electron Green's function Q e ^- U defined on 
the Matsubara discrete frequencies e n = nT(2n + I), 
uj n = 2irTn (T is the bath temperature) and expanded 
in series in the ac field A u . Any (local in time) 
observable is expressed through the sum TJ2 £ ^s,e-uj- 
The non-linear term of the fc-th order in this sum 
is E u b ^M 1) ,..^W)A kl ,i,...A tjW , where = 
lixT ri; and: 



(i) 



Tj2Go{ie)G (ie- 



)G (ie 



,(2) 



)... Go(is — iu>). 



(2.1) 



In the above equation we omitted the current opera- 
tors j or the position operators f coupled to the vector- 
potential A w or to the electric field E w = iwA w in the 
electron-field interaction: 



H e -f = 



-jA w transverse gauge 
-rE w longitudinal gauge 



(2.2) 



and introduced the exact electron Green's function in the 
absence of the time-dependent perturbation: 



G (ie) = }_^ , 



(2.3) 



where ty m (r) is an exact electron wavefunction in the 
mesoscopic system that corresponds to the stationary 
state with the energy e m . 

Causality requires the physical (retarded) response 
function K(u>x, ...Uk) which depends on the continu- 
ous frequencies Wj, to have no singularities in the up- 
per half-plain of each complex variable Wj. Thus 
K(uji, ...uik) is obtained by the analytical continuation 



) from the imaginary discrete points 
LJi = into the upper half-plane ^su)i > 0. In order 
to implement the analytical continuation one represents 



the sum over e n in Eq.(2.1) as a contour integral over the 
contour C that comprises all the points ie n = iirT(2n + l) 
(see Fig. la): 



x G (e)G (e - iaj w )G {e - iu: {1) - iuj {2) ) 



) (2.4) 



The next step is to deform the contur along the cuts 



0, 3e 



\5£ 



,(1) 



A 2 ' and so on (Fig. lb). 



In doing that one has to take care of the analytical 
properties of the electron Green's functions Gq{e): the 
integrand in Eq.(2.4) should be regular in each strip be- 
tween the neighbour cuts. This means that for 3e larger 
than that of the upmost cut all the Green's functions 
Go should be chosen retarded G$. In the strip just be- 
low this cut one Green's function with the argument that 
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takes zero value on this cut should be switched to the ad- 
vanced Gq. In the next stripe another retarded Green's 
function is switched from the retarded to the advanced 
one and so on. Thus for 3e < all the retarded Green's 
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be traced back to the causal nature of the non-linear re- 
sponse. 



B. The Keldysh diagrammatic technique 



FIG. 1. Contours of integration: (a) initial, (b) after de- 
formation along the cuts. The analytical properties [retarded 
(R) or advanced (A)] of Green's functions are also shown in 
each strip. 

The next step is to shift the integration from along the 
cuts to the real axis — oo < E < +00. To this end we shift 
the variable of integration e — » + corre- 



sponding to the (p + l)-th cut 



-J p \ Note 



that since iu>^ is a period of tanh(e/2T) the above shift 
of variables does not change tanh(e/2T) tax\h(E/2T). 
As a result we get in the r.h.s. of Eq.(2.4) a sum of the 
form: 



p=0 



E 
2T 

x G$(E - iuj^)...G^(E - i(co^ +1) + .../»))} . 

Now we are in a position to implement the analyti- 
cal continuation. It reduces simply to the replacement 
iu>W — > LUi because all the Green's functions have the an- 
alytical properties t hat g uarantee the regularity of each 
term in the sum Eq^^) in the upper half-plane of each 
of the complex variables Ui. 

One immediately notices the characteristic feature of 
Eq.(2.5): it is a sum of products (strings) of Green's func- 
tions such that each string is a product of p functions 
Gq followed by the product of k — p + 1 functions Gq. 
There is only one point in each string where the change 
of the character of analyticity occurs. This property can 



The analytical structure identical to Eq.( [2.5D arises in 
the Keldysh technique from the triangular matrix struc- 
ture of electron Green's functions Jll|[l4[ : 



G 



G R G K 
G A 



(2.6) 



where the superscript K stands for the Keldysh compo- 
nent that determines all the observables O: 



0(t) = iTr{dG K (i,t)} 



(2.7) 



Treating the applied ac field as a classical field, we assign 
Jl5| , [i4j the matrix vertex T°7i c -f to the electron-field in- 
teraction 7i c -f , where t° is the unit matrix in the 2x2 
Keldysh space. 

Using the usual expansion of electron Green's function 
G in powers of electron-field interaction: 



00 

0(t) = 1 Tr { 6 [{G H c _ f Y G ] K } (2. 

p=0 

i^^Tr[d{G^} l C%{1U- t Gt} p - 



P =o 1=0 
and the ansatz fi"^|~~i] : 



E \ 



+ + ...W p) ))...g£ceW p) ) (2.5) 

x [G$(E) - G$(E)] tanh 



G^E) = [G«(E) - G${E)] tanh j (2.9) 

one immediately obtains the same R-A structure as in 
Eq.©. 



III. NON-LINEAR RESPONSE IN THE TIME 
DOMAIN 



A. Non-locality of perturbation series in the time 
domain 



Let us consider the structure of the expansion Eq. 
TsJ) , for a given number p of the field vertices, in more 
detail. This will later help us to establish the structure of 
an essentially non-linear expessions for different observ- 
ables in the presence of a time-dependent ac field. 

Three different contributions can be distinguished 
there. In the time-domain representation, the first con- 
tribution can be depicted by the electron loop [Fig. 2(a)]. 
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(a) 



+ i / dr/(r) / df 



(b) 



(c) 



FIG. 2. Graphic representation of the expansion Eq. 
(|2.S|): (a) retarded, (b) advanced, and (c) retarded-advanced 
loops. Rays emanating from the electron loops correspond 
to the electron-field interaction TC c -{, and a triangle stands 
for the observable operator O. In the time domain, the fac- 
tor /(t) is assigned to the O-vertex in the diagrams (a) and 
(b). In the diagram (c) this factor is assigned to the re- 
tarded-advanced junction shown in bold. 



It entirely consists of retarded Green's functions: 



. dtpTr 



(3.1) 



O Gq (t, ti)H e -f(ti)G (ti, t2)H e -{(t2) ■ ■ ■ Go'(tp-i,t p )H e -f(tp)G^(tp, t — i 
The second contribution is associated with the electron loop [Fig. 2(b)] that solely contains advanced Green functions: 
-i I drfir) f dti ...d^Tr [dG^i + T,*!)^^ . (3.2) 



The third contribution is associated with the electron loop [Fig. 2(c)] built of I retarded and (p — I + 1) advanced 
Green functions: 



dt\ . . . dti-idti+i . . . dtpTr 



O Gq (t, ti)7i e -i(ti) . . . Go (i;-2, ^-i)W e -f (£z-i) 



drdr'fiT - t 1 ) G?(t,_i,T)[7<o_f(T) - W -f (t')]G£(t', t l+1 ) H c - f (t l+1 )G£(t l+u t l+2 ) . . .H e - f (t p )G£(t p ,t)] (3.3) 



where 



f(r) 



retarded— advanced junction 



27T smn(7rrJ ) 



oo 

G R ' A (r,r';t,t')=J2{[ G o' A ^ 



(3.4) 



p=0 



G 



R.A 



} , , 
J r. r';t,t' 



(3.6) 



denotes the Fourier transform of f(E) = tanh(_E/2T). 

The characteristic feature of the diagrammatic expan- 
sion for an observable in the time domain is that there 
is one special point (ray) on the loop of Fig. 2 which does 
not correspond to a single point in the time domain. It is 
the point where the Fourier transform of the energy dis- 
tribution function is assigned to. Of special importance 
is the retarded- advanced junction: 

G*(t, t) f(r - t>) [W e _ f (r) - We-f(r')] G A (r', t'). (3.5) 

It reveals the non-local in the time domain structure of 
the point [see Fig. 2c] where an (arbitrary long) sequence 
of retarded Green's functions is switched to an (arbitrary 
long) sequence of advanced Green's functions. 



B. Diffusons and cooperons in the time domain 

With the aim to describe the essentially non-linear de- 
pendence of observables on the external time-dependent 
field we introduce the infinite sequence of retarded (ad- 



where multiplication assumes the convolution over 
the coordinate and time variables {AB} r r >. t t > = 
J dt"dt" A r , r „ (t , t") B r „ >r - (*" , if) . 

In describing weak localization and mesoscopic phe- 
nomena a special role is played by the disorder averages 
of a pair of electron Green's functions called 'diffusons': 

V rr ,(t + ,t' + ;t-,t'_)=5(r ] -r ] ')D ri (t,t';r,r') (3.7) 
= (27n,r c )- 2 (G>, r'; t+, t' + )G A (v', r; t'_, *-)), 



and 'cooperons': 

C rr ,(t+,lf + ;t-,lf_) = y(t-t')C t (r), V ';r,r') 
= (27r^r )- 2 (G^(r, r'; t + ,t' )G A (r, r'; L,f_)). 



(3.8) 



vanced) Green's functions G 



R,A 



In Eqs.(3.7, 3.8) (...) stands for the disorder average, v 
is the mean electron density of states, r c is the electron 
momentum relaxation time, t± = t± n/2, t' ± = t' ± rj' /2. 
The 5-functions in Eqs.(3.7, 3.8) result from the assump- 
tion on the constant mean density of states over the 
relevant energy interval much smaller than the electron 
bandwidth. 

In the transverse gauge, the functions D v (t,t';r,r') 
and C t (r], 77'; r, r') obey the following equations jl6| which 
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correspond to the diffusion approximation with D being 
the diffusion constant: 



0_ 

at 



-y + D 



iV + A r it + 



V 



V 



electric field inside the system E(t) = —dA/dt is homo- 
geneous in space. 

In this situation the boundary conditions read: 



(3.9) 



S(r - r') 

xJ),(M';r,r')=^-f')4 T~- 

2-kvt?; 



and 



2 — + 7 + D 
on 



iV- 



xC 4 (77,77'; 



2S(i] — i] 



t-2 
2 



2-kvt?; 



where we assume the electron- field interaction Eq.( |2.2| ) 
corresponding to the transverse gauge with the weak 
space (on the scale of the elastic mean free path I — vpT c ) 
and time (on the scale of r e ) dependence of the external 
classical field A T (t) which is also supposed to be weak 
enough |A r (t)|.£ -C 1. We also assume the possibility 
for electrons to escape the mesoscopic system which is 
described by the (small) escape rate 7. 



C. The ring and the quantum dot geometry 



Equations (3.£, 3.10D should be supplemented by the 
boundary conditions. Below we consider two principally 
different geometries shown schematically in Fig. 3. 





FIG. 3. The ring (a) and the quantum dot (b) geometry 

One of them corresponds to the mesoscopic ring of 
the circumference L with a small aspect ratio pierced 
by the time-dependent magnetic field H(i) which gen- 
erates the circular electric field E(i) = —dA/dt, where 
A(t) = ±H x r. In this ring geometry the diffusons and 
cooperons should obey the periodic boundary conditions. 

The ring geometry corresponds to the pumping by the 
magnetic part of the microwave field when the size of 
the mesoscopic system L is less than the size of the skin 
layer. However, for a semiconductor quantum dot with 
not too large conductivity a the electric part of the mi- 
crowave field can also penetrate inside the mesoscopic 
system. For the usual situation L <C A where L is the 
size of the system and A is the microwave wavelength the 



iV + Q^M) D, 1 (t,t'-v,v') = 0, 

- II 

iV + Q^(t,n)} C t (r),T 1 ';v,v')=0, 



(3.11) 



where n is the vector normal to the boundary at a point 
r and we introduce a short-hand notation: 



(3.10) 



V 



n 



(3.12) 
(3.13) 



D. Diffusons and cooperons in the ring geometry 

In this case the convenient coordinate x along the ring 
is proportional to the azimuthal angle x = LBjlix. One 
can suppress the transverse coordinate, as the depen- 
dence on that coordinate at a small aspect ratio is negligi- 
ble. The tangential component of the field A x is indepen- 
dent on x. Given the periodic boundary conditions for 
diffusons and cooperons, one can switch to the Fourier- 



transforms D v (t,t';q) and C t (n,n';q) in Eqs.(p.9|, |3.10j ) 
with the quantized momentum q m — (2ir/L)m, where 
rn = 0,±1,±2.... Then the solution is straightforward: 



We"^-*') - D f tidT [ q - Q ^ ( ^)]\ (3U) 



D v {t,t';q) = - 



-f f> [,-Q<°> (t,r)] 2 



27WT? 



wher e t is th e step function and Q^ d,c ^ is defined by 
Eqs.([3T2| |3~i~5| ). 

An important particular case is the zero-mode diffu- 
sons D v (t, t') = D v (t, t';q = 0) and cooperons C t (n, n 1 ) — 
Ct(r], 77'; q = 0) which correspond to q = 0: 



t-t>e 
2'kvt'? 



D„{tJ) = - 
Ct{r),r]) = 



-D f dr [Q< d »(r,„)] 2 



(3.15) 



One can see that they decay as a function of t — t' or 
77 — rj' even at an escape rate 7 = 0. This is the man- 
ifestation of dephasing by the time-dependent external 
field. We note that Eq.(3.7) with n = corresponds 
to the electron density correlation function. In the ab- 
sence of electron escape the total number of particles is 
conserved and therefore Z3^ = o(^ t': 5 = 0) must be a con- 
stant for any t > t' . This is consistent with the property 
of QW(t, 77 = 0) = 0. However, there is no constraint 
that would prohibit a decay of D v (t, t') at a non-zero 77. 
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E. Diffusons and cooperons in the quantum dot 
geometry. 

This case is principally different because of the field- 



dependent boundary conditions Eqs.(3.11) and the po- 
tential (longitudinal) nature of the electric field inside 
the dot. This makes the description in the longitudinal 
gauge more convenient in the quantum dot geometry. 



1. Equations for the diffusons and the cooperons in the 
longitudinal gauge 

Performing the gauge transformation: 
D v (t, f; r, r') = e*^" D v (t, t'; r, r') e -iQ w V,vW t 



(3.16) 



<. r ^ = P ,Q (c »(M)r£ t fa v , ;r)r ') e -iQ M My 



C t (ri,ri ;r,r ) = e 



(3.17) 



we switch to the longitudinal gauge. 

This transformation removes the time-dependent field 



from the boundary conditions Eq.(3.11). However, 
it remains in the equations for D v (t,t';r,r') and 
Ct(?7, 77'; r, r') but only as the time-derivatives: 

Q(t,v) = ^Q (d) (i^) = 2^-Q^(t, v ) (3.18) 

-.(,-§)-.(J§), 

where E(t) = -dA(t)/dt. 

The operators in the l.h.s. of the corresponding equa- 
tions take the form: 



C d = \—+ 1 -DV 2 + irQ(t, V ) 



drj 



7 - DV 2 + irQ(t,r]) 



(3.19) 



(3.20) 



It is convenient to expand D v (t, t'; r, r') and 
C t (r),Tj';T,r') in an infinite sum A Vfl ^(r)^^') 

over eigenfunctions § v (r) of the diffusion operator — Z?V 2 
with the Neumann boundary condition [Eq.(3.11) with 
q(<2,c) _ gj ^hen ^ ne equation for the amplitudes A vyL 
is of the form (for simplicity we consider the case 7 = 
and suppress the redundant indices and variables): 



d_ 

&t 



A Ufl + E u A vit + iQ(t) r ^ A ^ = *(* - 0* 



(3.21) 



where r^ is the matrix element of the vector r in the 
basis of $„(r) and E v is the eigenvalue that corresponds 
to the eigenfunction $„(r). 

Though this equation looks similar to the (imaginary 
time) Schrodinger equation for the dynamical Stark ef- 
fect, identifying Q as the actual electric field inside the 
dot would lead to a mistake. Let us consider an impor- 
tant case of the constant in time electric field E which 
corresponds to A(t) = -Et. Then Eqs.Q3.12j, [Tl|) give 
Q( rf )(t,r?) = - E77 and Q( c ~>(t,r)) = -2Et. One can see 
from Eq.(3.18) that in both cases the corresponding time- 
derivatives Q( d ) and Q' c ' are identically zero! 

The conclusion which can be immediately drawn from 
this observation is that the constant in time longitudi- 
nal (potential) electric field cannot lead to a dephasing 
[cf. Rcf. fril l. This statement is not true for the circu- 
lar electric field considered above. This field is not con- 
stant in the Cartesian coordinates and is not potential 
curl E^O. 



2. The weak-field adiabatic and anti-adiabatic limits 



The general solution to Eq.(3.21)is unknown even for 
the space- homogeneous electric field E(i) and in the er- 
godic limit 



E c (t-i')>l E c {f]-rf)>l. 



(3.22) 



However, one can find a simple approximate solutions in 
the weak field limit 



\Q\L 

E c 



< 1. 



(3.23) 



In the ergodic limit Eq.( |3.22 ) the diffusons D v (t, t'; r, r') 
and cooperons Ct(r), if; r, r') are nearly space- 
independent, i.e. the corresponding expansions are dom- 
inated by the it zero-mode $0 with Eq = 0. By definition 
the next mode has the eigenvalue equal to the Thouless 
energy E\ — E c . Fo r the corresponding amplitude Aqq 
we get from Eq. ( 3.21 ) : 



d 

—A 00 + iQ(t) 22 r Qfl A^o = S(t - t'), 



(3.24) 



where we choose the system of coordinates in which 
roo = 0. 

In the weak field limit Eq.( |3.23 ) one can neglect in 
Eq.( |3.21 ) Qr^AyQ compared to E^A^o {l*,v 7^ 0) and 
obtain the closed system of equations for Aqq and A^q: 



"a 

dt 



EpAfti + «Q(i)i> A 00 = 0. 



(3.25) 



Solving Eq.( 3.25| ) and substituting A^q into Eq.( 3.24 ) we 
obtain: 
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9 Ann + f Q(t)A(* - t")Q(t") A m {t",t') dt" = S(t - t'), 



dt 



(3.26) 



where A(i) is the matrix in the vector space: 

[AWl^EW-^'ia- (3-27) 



the (5-function [A(i)]y = CijS(t), where: 



(3.28) 



Eq.( |3.26 ) can be further simplified in the adiabatic limit 

where Q(i) is dominated by th e frequencies u> <C E c . In Then Eq.( 3.26| ) can be immediately solved and we obtain 



this case the quantity Eg. ( 3.27 ) can be approximated by for the zero- mode amplitudes: 



A$(t,t';r ] )=D ri (t,t') = 



1-kvtI: 



exp 



(3.29) 



A$(r ] ,r, , ;t) = C t (r 1 ,T 1 ') = 



"tj—tj' 

IuvtI: 



exp 



in Eqs.(3.15) contains a quadratic form in Q and is in- 



dependent of L. 



F. Time-dependent random matrix theory (TRMT) 



QfWOaiQf (t,v")dr, 



(<=), 



Eqs. (F29|) ,(|3~30|) are valid in the adiabatic limit lu <C 
E c provided that the conditions Eqs.( 3.22 ),( 3.23 ) are ful- 
filled. 

In the opposite anti- adiabatic limit lu E c , the in- 
tegral term in Eq.( 3.26 ) is recast as the total time- 
derivative (which has zero time average and thus strongly 
oscillates) and the remainder consisting of the strongly 
oscillating term — Q(i) A.(0)Q(t)A oo (t, t') and the term 



dt" Q(i) [d/dtA(t-t")]Q{t")A 00 (t",t'). 



that contains a weakly oscillating part. Integrating 
this term by parts and using the identity A(0) = 
d / dt[A]ij\ t= o — ~DSij one extracts this weakly oscil- 
lating part: 

- Q(t) A(Q) Q(t) Ao (t, t') = D [Q(t)} 2 A 00 (t, t'). (3.31) 

Not surprisingly, we conclude that in the anti-adiabatic 
weak-field limit the system does not feel the boundary 
and the zero-mode diffusons and cooperons in the quan- 
tum dot geometry coincide with those in the ring geom- 
etry Eqs.(|3~15|). 

However, in the adiabatic weak-field limit these two 
geometries are principally different, since the dephasing 



factor in Eq.(3.2E),(3.30) contains a quadratic form in 



the time-derivative Q and a structural constant that de- 
pends on the system size L, while the dephasing factor 



The zero-mode approximation Eq.( 3.22[ ) is equivalent 
to the random- matrix theory (RMT). The fact that there 



'/' 



(3.30) 



are two different forms, Eqs. (|3.29| ),(^30|) and Eqs. (|3.15|) , 
of diffusons and cooperons in the zero-mode approxima- 
tion suggests two different ways of defining the time- 
dependent RMT. The idea is to define the Gaussian en- 
sembles of time- dependent random matrices which repro- 
duce the expressions Eq. (|3~29l) , ( |3~30| ) or Eqs . (flj) . Then 
so defined TRMT can be applied to describe not only dis- 
odered mesoscopic systems with the diffusion motion of 
electrons (used in the above derivation) but also ballistic 
quantum dots with the chaotic electron motion. 

Before proceeding with the formal derivation we ad- 
dress a possible confusion based on the common wisdom 
that only systems in the external field with the charac- 
teristic frequency lu <C E c can be described by the ran- 
dom matrix theory. This statement is valid for a linear 
response but it is incorrect in the non-linear case. The 
point is that in the linear case the frequency of the exter- 
nal field enters all diffusons or cooperons as [Dq 2 — ito]" 1 
where lu is the difference between the energy variables of 
retarded and advanced electron Green's functions. As- 
suming the summation over momenta q and the fact that 
the first non-zero mode corresponds to Dq 2 = E c one 
concludes that at u> <C E c the main contribution to the 
sum over momenta is given by the zero mode with q = 0. 
In the non-linear case the situation is more complicated, 
since TL c -{(t) in Eq.( |3.6D is a sum of two parts propor- 
tional to e luJt and e~ lut . As a result of the frequency fu- 
sion lu — lu = in the field-dependent diffuson(cooperon) 
self-energy part the difference between the energy vari- 
ables of retarded and advanced electron Green's func- 
tions constituting a difuson(cooperon) may be zero de- 
spite lu 3> E c . This is the reason why the high-frequency 
external field modifies the zero-mode approximation but 
does not kill it. The two modifications of the TRMT are: 
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(i) For the conductors of the ring topology and the 
quantum dots in the anti-adiabatic limit, uj 3> E c , the 
time dependent RMT is defined by the matrix Hamilto- 



/ Tjnra rrn m 
\ H H 



NS 2 

7T 2 

sc, 



H = H + i£(t)V a , 



(3.32) 



where £(4) is a real-time dependent function, N x N ran- 
dom matrix Hq belongs to the Gaussian orthogonal en- 
semble while V a is the real anti- symmetric random ma- 
trix (which corresponds to the sign — in Eq.( 3.33| )) of the 



(K7C m ) = -[^'^' T <W<W], (3.33) 
7T 



(5 being the mean level spacing, and C is a non-universal 
constant to be identified later on. To make a link be- 



same size, 



tween the TRMT, Eqs. fl3.32p and ( |3.33| ), and the zero- 
dimensional approximation Eqs. (3.15) and (3.31), it is 
useful to introduce the TRMT diffuson, D n (t,t'), and 
cooperon, 64(77,77'), propagators 



^2{G% m (t+,t' + )Gi n (t'_,t-)) = (27nyT c (N)) 2 S(n-r/)D TI (t,t'), 
J2( G *™(t+A)G* m (t^l/_)) = l -(2irvT c {N)f5(t-t')C t {n 1 n'), 



(3.34) 



where G R > A (t,t') are retarded or advanced Green's functions that correspond to the matrix Hamiltonian H, t± 
t ± 77/2, 4 = t' ± n'/2; t c (N) = ir/(2N6), and v = 1/6. Using the standard method of Refs. |l6| (see also Ref. | 
we derive the following equations for the TRMT propagators: 



d 



C 



dt 

d_ c r 

d^ + ~2 



S(t-t') 



f( t+ !)-«(* -I)] h^-s^w 
{ ( i+i ) +f ( t -l)]n c ,,„',-*-') 



2ttvt%(N)' 



(3.35) 
(3.36) 



Comparison with the microscopic Eqs. (3.9) and (3.10) 
[or (3.31)] suggests identifyi ng th e phenomenological con- 
stant C introduced in Eq. ( 3.33| ) with the diffusion coef- 
ficient D; the time-dependent function £(4) plays the role 
of the vector potential A(t). This establishes the equiva- 
lence between the TRMT and the zero-dimensional limit 
of the microscopic theory of Sec. Ill on the perturbative 
level. 

(ii) For the quantum dots in the adiabatic limit, uj <C 
E Cl the time dependent RMT is defined by the matrix 
Hamiltonian p|,|ll[ 



H = H +C(t)V s 



(3.37) 



where £(4) is a real-time dependent function, H and V s 
are statistically independent, N x N, random matrices 
belon ging to the Gaussian orthogonal ensemble [ see Eq. 
( |3.33 ) with the sign +]. 

Ag ain, a link between the TRMT, Eqs. fl3.37p and 
( |3.33| ), and the zero-dimensional approximation Eqs. 
(3.29) and (3.33) of the full microscopic theory, is easily 
established by deriving the equations for TRMT propa- 
gators. One obtains [see also Refs. 



d_ 



C 
2 



(5(4 - Q 
27ri/r e 2 (^)' 

gfo? - v') 

2ttvt2(N)' 



(3.38) 
(3.39) 



Comparing with the microscopic Eqs. (3.29) and 
(3.30), we co nclude that the non-universal constant C 
of Eq. ( HI) has to be identified with the one given by 
Eq. (3.28); the function £(4) plays the role of ac electric 
field E(t) in Eq. (3.18). Thus, equivalence between the 
microscopi c app roach of Sec. Ill and the TRMT of the 
formEq. (^37p is proven perturbatively. 



IV. THE 'LOOSE' DIFFUSONS 

In this section we show that starting from the 
quadratic in the external field order the diagrams of the 
impurity technique acquire a new feature: one can draw 
the diffuson with the free end ('the loose diffuson') which 
carries zero momentum and zero frequency. We will show 
below that it is exactly the element which describes heat- 
ing by the external field. 
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A. The loose diffusons and the retarded-advanced 
junctions 



1. Loose diffusons in the ring geometry 



The analytical structure of the retarded-advanced 
junction Eq^JiJj) leads after the disorder averaging to 
an unusual object, the 'loose' diffuson. Let us consider 
this object for the simplest ring geometry. 



Consider the part of a diagram for an observable or 
a product of observables that contains the retarded- 
advanced junction (Fig. 4a). 



R A A 



R R A 



R A 




(a) 



(b) 



(c) 



(d) 



FIG. 4. The retarded-advanced junction (a), and the Green functions' disorder averaging [(b) and (c)] resulting in the 'loose' 
diffuson shown by the wavy line (d). It ends up with the triangle; one of the triangle apexes is built upon the retarded-advanced 
junction, see (b) and (c). The diagram (d) corresponds to the disorder averaging depicted in (b). Dashed lines in (b) and (c) 
denote the diffuson. Because of the vector nature of the vertex TC e -f in the transverse gauge Eq.(2.2) the loose diffuson should 
'embrace' two vertices 7G-f in order for the integral over the directions of electron velocity to be non-zero. 



One can isolate the retarded-advanced junction from the 
rest of the diagram by performing the disorder averaging 
as shown in Fig. 4b, c. As the result the 'loose diffuson' is 
formed (Fig.4d) which originates from the main body of 
the diagram and terminates at a triangle that consists of 
the retarded-advanced junction and another field vertex 
adjacent to it. Given the condition |A r (i)|.£ 1 one can 
neglect the loose diffusons terminating by a polygon with 
the number of field vertices larger than two. 

We stress that the loose diffuson can be built only 
using the retarded-advanced junction. Indeed, the tri- 
angle in Fig.4d whose edges correspond to the aver- 
age retarded and advanced Green functions G^ A ^ (p) = 



Jdr(G^ A) (r,r';E))e-^-^ = (E - £ p ± i/(2r e ))- 
describes the electron motion on the ballistic scale with 
the electron momentum relaxation time r being the 
smallest time scale in the problem. One may effectively 
approximate the average Green's functions in the time- 
momentum representation G R( - A ' (p, t) by the <5-functions 

of the form G R < A )(p,i) ~ G^> (p)S(t). With this ap- 



proximation and Eq.(3.5) the triangle in Fig.4d reduces 
to the quadratic in A(i) combination 



[A^-A^Aft) /ft-*!) 



multiplied by a constant (d is the dimensionality of mo- 
mentum space) 



(v%/d)v / d£(p)Gl 



=o(P) [Gk 



, (p)} 2 = 2mvDr e . (4.2) 



The triangle that corresponds to Fig. 4c is given by: 

- 2mvDr e [A(t[) - A(*i)] A(«i) f(t[ -h). (4.3) 

Without the retarded-advanced junction, all the 
Green's functions in the triangle would have the same 
analyticity (R or A), and the integral over £(p) vanishes. 

Since in the dynamical approach (e-e and e-ph inter- 
actions are neglected) there is at most on e re tarded- 
advanced junction per electron loop [see Eq.( |3.3| )l, there 
could be not more than one loose diffuson for a dia- 
gram describing the disorder-average of a single observ- 
able (0(t)) and not more than k loose diffusons for a 
diagram describing the disorder-average of a product of 
k observables (0i(ii)...O/s(ifc))- 

Each of them is given by: 

T>(t, rj) = 2mvr e D J dt x J drx[A{tx + rj) - A(ii)] 2 

x/^ZL^ + ^ti + ^n-r), (4.4) 



(4.1) where -q = t' — t = t[ — t\ 



and both contributions 
Eqs.(4.2),(4.3) to the triangle have been summed up to- 
gether. In the ring geometry only the zero mode part of 
D v (t, t';ri — r) survives integration over the coordinate 
rj of the free end and using Eq.(3.15) we finally obtain 
V(t,ri)=iT- 1 f(ri)A r ,(t), where: 



A ri (t) = D J d£[A(£ + ??)-A(£)] 2 e-^-«) expj-Dy d£'[A(£' + rj) — A(£')] 2 j. 



(4.5) 
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2. Loose diffusons in the quantum dot geometry 



as follows: 



Eq.(4.5) holds also in the quantum dot geometry in 
the anti-adiabatic case u) 3> E c . In the opposite adiabatic 
case u) <C E c the longitudinal gauge is more convenient 
than the transverse gauge originally accepted in the pa- 
per. Eqs.(^29|),(^30|) are nothing but the diffusons and 
cooperons in the longitudinal gauge in the zero-mode ap- 
proximation. Using Eq.(3.5) with 7i c _f corresponding to 
the longitudinal gauge Eq.(2.2) the loose diffuson can be 



represented in terms of the matrix elements Affl and r 



2?(t, rf) = 2mvr*f{ri) / dt x [dA{t x + rf)/8t x - dA(t x )/dt x ] 

T) 



i 1- 2 ; ^ 7] l *>' 



(4.6) 



where A$ obeys Eq.(|Uf) with Q(t) = Q(t, -rf) and 



A$ given by Eq.(|3.29|) 



Solving this equation and using the (^-function approxi- 
mation for A(f) we obtain T>(t, rf) = ir^ 1 f(r))A v (t) with: 



A„(t) 



dZC\E(t + T))--E(£)] 



2 e -r<*-« exp 



decme+v)-na? 



(4.7) 



where C[E(£ + rf) — E(£)] 2 is the short-hand notation for 
[E(e + r y )-E(0] l C iJ [E(C + ? 7 )-E(e)] J and E = -dA/dt 
is the time-dependent electric field. We also re-installed 
the finite escape rate 7. 



B. Loose diffusons and the singularity of the 
quadratic response 

Eqs.(|4.5|).( |4.7| ) have similar structure: A v (t) contains 
a quadratic in the external field pre-factor multiplied 
by the exponential dephasing factor. Because of the 
quadratic in A pre-factor the loose diffuson does not arise 
in the linear response theory. However, if one considers 
the quadratic in A response, the loose diffusons must be 
taken into account while the field-dependent dephasing 
should be neglected. In this approximation we have: 



A„(i) 



d^[A(e + r;)-A(0] 2 e- 



-7(*"«) 



(4.8) 



and the similar expression in the quantum dot geometry. 

Consider the harmonic pumping A = A cos(wi) 9{t) 
that is switched on at t = 0. Then at a time t ^> co^ 1 
the loose diffuson averaged over the period is given by: 



A v (t) 



2DAl 



1 



--ft 



7 



(4.9) 



One can see that in the limit 7^0 the loose diffuson is 
linear in t. 

We have already mentioned that the disorder average 
of the product of k observables contains at most k loose 
diffusons. Were the linear in t growth in Eq.(4.£) un- 
restricted, this would mean that the typical value of a 
mesoscopic observable grows linearly with the running 
time. For a particular case of the direct current arising 
in a mesoscopic ring under ac pumping this statement 
can be found in Ref. pM. 



Another similar statement concerns the ste ady-state 
regime when the limit t — > 00 is taken in Eq.(4.9) prior to 
the limit 7^0. One could argue from Eq.(4.9) that the 
typical value of a mesoscopic observable in the steady- 
state is diverging as 7 — > ! 

However, Eqs . (|4 . 5|) , ( |4 . 7|) clearly show that both state- 
ments are artefacts of the quadratic approximation. In 
fact because of the field-induced dephasing the quantity 
A v (t) defined by Eqs. (4. 5), (1. 7) is always smaller than 1. 
What is really singular in the limit 7 — > is the quadratic 
response susceptibility. However, it does not mean a di- 
verging mesoscopic quantity, since at 7 — > the region of 
validity of the quadratic in A approximation shrinks to 
zero. 



C. Loose diffusons and the electron energy 
distribution 



Note that originating from the retarded-advanced 
junction Eq.( |3.5| ) the loose diffuson is proportional to 
the combination f(rj)A ri (t), where f(rj) is essentially 
the Fourier-transform of the Fermi distribution function. 
From the procedure of building the loose diffuson it is 
clear that any diagram with the loose diffuson has a 
parent-diagram without the loose diffuson. For a partic- 
ular case of the variance of persistent current in a meso- 
scopic ring the diagrams of Fig. 7 (or Fig. 8) with one or 
two loose diffusons stem from the diagram of Fig. 6a (or 
Fig. 6b) that contains no loose diffusons. One can check 
(see Sec.V) that the sum of all diagrams of the given fam- 
ily is equivalent to replacing f{rj) in the parent diagram 
by: 



/ , (t + 7?,i)-[l-A,(i)]/(r/). 



(4.10) 



Consider now the simplest case of the steady-state. It 
corresponds to the ac pumping switched on at t = -00. 
Let us define the time average: 
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D[A(t + rj) - A(t)} 2 ring or dot with uj > E c 



C\E(t + v)-E®] 2 



dot with oj <C E c 

(4.11) 



Substituting the time average T(rj) for D[A(t+rj)— A(t)] 2 
or C[E(t + v)~ E (*)] 2 in Eqs.(g|),(|J) we obtain the 
function F(t + rj, t) = F(r]) that depends only on the 
difference of its arguments: 



F(V) = 



7 + r(77)' 



(4.12) 



We conclude that the loose diffusons amount to the 
renormalization of the energy distribution function in a 
parent diagram. 

This observation is an example of a generic rule that 
the 'loose propagators' in a filed theory can be elimi- 
nated by the proper choice of the initial state ('vacuum') 
that in kinetics includes also the energy distribution func- 
tion. Eq.([yJ) gives the form of this distribution for an 
open dynamical system with no intrinsic relaxation and 
with the electron escape rate 7. For the particular case 
of a harmonic pumping in the quantum dot geometry 
Eq. ( [4.12D has been established in Ref . Q . 

The renormalized energy distribution function 



Eq.(4.12) retains the property F(E = ±00) = ±1 
of the equilibrium distribution f{E) = t& nh(E /2T). 
This follows from the fact obvious from Eq.( [4.1l| ) that 
T(rj — > 0) = 0. However, the form of the energy distribu- 
tion is different from tanh(i?/2T) and strongly depends 
on the spectral content of the pumping field. For low bath 
temperature T it contains at least two energy scales: the 
bath temperature T and an additional scale set by 
the condition r(l/T*) = 7. For a harmonic pumping 
r(?y) = 47/V sin 2 (w7//2), where: 



DE(t) 2 



CE(t) 2 

7 



ring or dot with to E c 
dot with lu <C E c 



For a strong pumping with M 3> 1 one finds 



(4.13) 



(4.14) 



This result |Tl[] corresponds to the diffusion in the en- 
ergy space with T» being the displacement and Af being 
the number of random- walk steps for the time 7 -1 each 
with emitting or absorbing the energy lu. In this case the 
inverse Fourier transform F(E) of F(rj) is dominated by 
the time intervals near zeros of T(rj) and we obtain: 



F(E) = 



2VN 



+ 00 

£ 

k— — 00 



tanh 



E 



2T 



(4.15) 



At T < lu < the function F{E) changes from -1 to 
+1 over the scale T* which plays a role of the effective 



electron temperature. However it has a fine structure of 
sharp small steps with the width of the transition regions 
being equal to the bath temperature T and the width of 
the plateaus being equal to to (see Fig. 5). 




FIG. 5. The anti-symmetric part F(E) = 1 - 2n(E) of 
the electron energy distribution n(E) for an open mesoscopic 
system of free electrons under harmonic pumping. Solid line 
corresponds to TV = 25, u) = 1, T = 0.1; the dotted line is 
tanh(£/2T»). 

In case of the white-noise pumping T(rj) = 27V7 = 
const for all rj ^ and we formally obtain — > 00. The 
same result follows from Eqs.( 1.14 ),( [4.13| ) in the limit of 
a closed system 7^0. On the diagrammatic level this 
manifests itself in the identity = 1 that holds in the 
limit 7 — > for all rj 7^ and causes a cancellation of all 
the diagrams in the given family ]l9| . 

In the absence of dissipation or an electron escape 
the result that the effective electron temperature in the 
steady state T* = 00 is trivially correct but is certainly 
unphysical. It clearly sets the limit of the dynamical 
approach and shows a necessity to take account of an 
intrinsic dissipation in closed mesoscopic systems. 



D. Loose diffusons and the kinetic equation 

In this subsection we demonstrate that the function 
F{t\,t2) defined by Eq.(llC) is indeed the solution to 



the kinetic equation. For simplicity we consider the case 
of the ring geometry and set 7 = 0. 

To derive the kinetic equation we start with the left 
and the right Dyson equations for the 2x2 Greens func- 
tion G: 



(G ( 7 1 -S)®G = <5(x 1 -x' 1 ), 
G^(Go 1 -S) = (5(x 1 -x' 1 ). 



(4.16a) 
(4.16b) 



Here, the standard notations 
x fc = (r k ,t k ), and 



L4[ were adopted with 
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G X ( x ) = i -Q t ~ £p=-iV r -A(t)- 



(4.17) 



As the external AC field is a classical field, the com- 
ponents £ K ( R ' A ) of the self energy 



E(r*,r'*' 



E A 



(4.18) 



e k ( r,a) = (j A )(x 1 )G K ( R - A )(x 1 ,x 2 )(]A)(x 2 ). (4.19) 



Subtracting the two Dyson equations Eq. (4.16) from 
one another, and taking the Keldysh component of the 
result, one obtains: 



G - 1 (x 1 )-G - 1 *(x 2 )l G K (x l7 x 2 ) = [S K ® G A — G R ® S K + £ R ® G K — G K ® S A ] (x 1; x 2 ). 



(4.20) 



Before doing the next step we introduce the time- 
dependent energy distribution function n(E,t) = 
(1/2) [1 - f(E, *)], where f(E, t) is related to the Keldysh 
component of the matrix Green function: 



G K (r 1 t 1 ,r 2 t 2 ) = j dT[F(t 1 ,T)G R (r 1 T,r 2 t 2 ) 



G A (r 1 t 1) r 2 r)F(r,t 2 



mt) = J ge"^( 



2' 2 



(4.21) 
(4.22) 



Now we substitute Eqs. (4.19) and ( 4.21 ) into Eq. 
(4.20), and perform the disorder averaging. Using the 
identity: 



F(t ll f 2 ) = ^ ; ^(G K (p;i 1 ,t 2 )) ! 



(4.23) 



we arrive at the equation 
d d 



0*1 5*2 



F(*l,*2 



-D [A(*i) — A(* 2 )] F(*i, * 2 ). 

(4.24) 



The latter is easy to solve by introducing the function 
F(t,rj) = F(t + n/2,t - 77/2), with * = (* x + * 2 )/2 and 
77 = *i — * 2 being the Wigner variables: 



A(« + ?)-A( t -| 



F(*,77). 
(4.25) 



In accordance with Eq. ( 4.22 ), the variable t has a clear 
meaning of the global running time. We supplement 
Eq.(4.25) by the initial condition F(t = —00,77) = fiv) 
to end up with: 



F(t + V ,t) = f( V )exp\-D 



d£[A(£ + »7)-A(0] 

(4.26) 



Comparison of Eq. (4.26) with (4.5) (at 7 = 0) shows 
that the function F(t±, * 2 ) is identical to the function 
F(t±, * 2 ) defined by Eq.(4.1C) which follows from the eval- 
uation of diagrams with loose diffusons. 



This proves that the loose diffuson A IJ (*) determines 
the time- dependent electron energy distribution function 
n(E,t) = (1/2) [!-/(£?,*)]: 



^e-^/(77)[l-A,(*)], (4.27) 

Z7T 



where A^(*) is given by Eqs. (4. 5), (4. 7) 



V. FLUCTUATIONS OF PERSISTENT 
CURRENT IN MESOSCOPIC METALLIC RINGS 
IN AND OUT OF EQUILIBRIUM 

For illustration purposes we consider in this section 
how the general formalism described above works in the 
particular problem of the persistent current |20[| in meso- 
scopic rings pierced by a constant magnetic flux <f> and 
subject to ac pumping. Since the disorder-averaged per- 
sistent current of non-interacting electrons considered in 
the grand-canonical ensemble is exponentially small H] 
we concentrate on the mesoscopic fluctuations of persis- 
tent current at temperatures T» 1 /tj, where 1 /r^ is the 
total dephasing rate including that of the ac pumping. 
This condition allows to neglect dephasing everywhere 
but in the loose diffusons which describe the evolution of 
the electron energy distribution under the ac pumping. 
We will show by straightforward diagrammatic calcula- 
tions that Eq.(4.10) indeed holds if all diagrams of the 
given family are taken into account. This illustrates how 
the diagrammatic technique takes care of the correct elec- 
tron energy distribution in the non-eqilibrium problem. 



A. Equilibrium fluctuations of persistent currents 

We start with the equilibrium fluctuations of persis- 
tent currents in order to specify the parent diagrams for 
the problem considered. 

Following a standard route, we express the persistent 
current in terms of exact retarded and advanced electron 



Green's functions G 



R(A). 



'PC 



Ug$(e)-g£(e)) 



(5.1) 
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Then the variance of the persistent current fluctuations 
is given by: 



(I: 



pc) 



= 2 



dE 



dE' 

~2T 



f(E)f(E') 



(5.2) 



To perform the impurity averaging in Eq. (5.2), it is con- 
venient to use a representation |2^j2j| in which slow diffu- 
sion and fast ballistic modes are explicitly separated from 
each other. The lowest order (single-loop) diagrams con- 
tributing to the persistent current fluctuations are shown 
in Fig. 6. There, wavy lines correspond to Diffuson 
(T>k) or Cooperon (Ck) propagators, while triangles and 
a square (whose edges correspond to the average Green 
functions) represent the electron motion on the ballistic 
scale. On the diffusion scale, the latter are reduced to 



certain constants. 




(b) 



FIG. 6. Local (a) and nonlocal (b) single-loop diagrams 
for the equilibrium persistent current fluctuations. Triangu- 
lar vertices mark the current operators j a . These are the 
parent diagrams for the persistent current fluctuations under 
ac pumping. 



Applying conventional rules (e.g., see Ref. p4|), 
convert the diagrams in Figs. 6(a) and 6(b) to 



we 



(i? c )i = ^vDrle 2 f ^/ 2 (??)^[C fem _ 2 ^ o (r ? ) -V km (r,)}, 

J m 

(I 2 c ) ni = (4nvDr 2 ) 2 e 2 f d V f 2 (rj) f d£ ]T [k 2 n V km {£)V km (r, - £) - k 



(5.3a) 



m-2d>/<t> C k m - 2 4, Un (QCk m _ 2Mn {v - £)]• ( 5 - 3b ) 



Here k m = (2n/L)m with m running over all integers 
represents the spectrum of diffusion modes allowed for 
conductor with the ring topology, L being the circum- 
ference. In Eqs. ( |5.3|) , the contributions of single- and 
double-diffuson (cooperon) diagrams have been singled 
out in the time domain, in which 



V k {t) = Ck{t) 



cxp(-Dk 2 t). 



(5.4) 



In fact, Eq. (5T) allows us to effectively compact- 
ify the nonlocal diagrams, Fig. 1(b), so that the total 
fluctuations (/p C ) are solely expressed in terms of the 
contributions of the local diagrams, Fig. 1(a): 



(4 C ) =8^Dr 2 e 2 J d V f 2 ( V ) [ D 



d 



d 



dD dD 
Vk m (v)}- 



D 



(5.5) 



The fluctuations of persistent currents, Eq. (5J3) are 
manifestly periodic in the flux 0, with the period 4>q/2. 
This can explicitly be displayed by performing the re- 
summation in Eq. ( |5.5| ) using the Poisson formula: 



n=l 



with 

{il) 



4C 2 n 2 I e 



T 



dx 



-n 2 /{ix) 



c 3 / 2 sinh 2 (7rfa;)' 



(5.6) 



(5.7) 



where T = T/E c denotes the electron temperature mea- 
sured in the units of the Thouless energy E c = I/tq. For 
generality we also introduce the coefficient C which is 
equal to 1 for the case of potential disorder (orthogonal 
ensemble) considered here and C = 1/2 for the strong 
spin-orbit interaction (symplectic ensemble). 



Eq.(5.7) is in a complete correspondence with the ear- 
lier results obtained in Refs. fMEfl. 



B. Effect of ac pumping on the persistent current 
fluctuations 



Now let us assume that a time-dependent circular 
field A(t) [see Fig. 3a] is applied to the mesoscopic ring 
and consider the persistent current fluctuations under ac 
pumping. The dc current in a ring (overline means the 
time averaging) Id c 



JW (i) + J(2) (i) can be found using 
Eqs.(|24-(M) [see also Fig. 2]: 



■ J dr, f(rj) Tr {jG*(f, t-r,)- iG A (t + V , t) } 



(5.8) 



I {2) {t) = ijdt 1 J d v f( v ) [A a (h+r}) - A a (h)} 

x TrjjG^i, ii + r])j a G A (t 1 ,t)^ , (5.9) 
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where G R ' A are exact electron Green's functions in the 



persence of ac pumping Eq.(3.£). 

One can check that for equilibrium electron Green's 



functions Eq.(5.£) reduces to Eq.(5.1). The contribution 
Eq.(5.9) is present only under ac pumping. It describes 
two principally different effects of ac pumping. One is the 
rectification of ac field discussed in Rcfs. |9||(|. This effect 
is similar to the photovoltaic effect in a single-connected 
geometry Another one is the heating by ac field 

which wc will study starting from the simplest zero order 
in the pump field parent diagrams of Fig. 6. The heat- 
ing effect in rectification (or photovoltaic effect) can be 



studied in a similar way starting from the parent di- 
agrams of the second order in the pump field pQ,0,E^. 



The daughter-diagrams with loose diffusons which 
arise after disorder averaging and correspond to the 
parent-diagram of Fig. 6a. are given in Fig. 7. We stress 
that although the parent diagram of Fig. 6a arises as the 
result of the disorder averaging of (if), the daughter- 
diagrams always involve I2 that contains the retarded- 
advanced junction and allows to build the loose diffuson. 
In a similar way we get the non-local daughter-diagrams 
that correspond to the parent diagram of Fig. 6b. 




(a) (b) (c) 



FIG. 7. Local single-loop daughter-diagrams with one (a) or two (b) loose diffusons contributing to the non-equilibrium 
fluctuations of persistent current. The diagram (a) arises from the disorder average (Ii(t) 12(f)) and the diagram (b) arises from 
the disorder average (I^ifjl^it')). They are complementary to the parent-diagram of Fig. 6(a). Hikami box of the diagram (a) 
is detailed in (c). 





□ = □ + □ 



FIG. 8. Nonlocal single-loop daughter-diagrams: the diagrams (a) and (b) arise from the disorder average of (Ii(t)l2(t')) 
while the diagram (c) arises from the disorder average {I2{t)l2(t')} . These are complementary to the parent-diagram of Fig. 
6(b). The corresponding Hikami box (d) is also shown. 



Calculating the diagrams of Fig. 7 and Fig. 8 and as- 
suming T»7, we neglect both dephasing and elec- 
tron escape in the loop diffusons/cooperons and adopt 
Eq.(5.4) to describe them. Then summing up all the di- 
agrams of Fig. 6-8 we arrive at the expression for the dis- 
order average (Ipc(t)Ipc(t')) which is exactly the same 
as Eq.(|J) with f 2 (rj) replaced by f 2 (rj) (1 - K n (t)){\ 



A- n (t')). This leads to the ansatz Eq.([4T0|). 

In particular, for a harmonic pumping with the fre- 
quency lo = ujE c we obtain from Eq.(4.12): 



AC 2 n z 
nV 2 



where 



F(x) 



TD 



T 



dx 



smh(nfx) [l + 4JV sin 2 (^)] 



-n 2 /(4x) F 2^ ( 51() ) 



(5.11) 



Here we assume that the ring is connected to the elec- 
tron reservoir by a passive lead [see Fig. 10b] which results 
in a finite electron escape rate 7 and allows to reach a 
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steady-state regime. The escape rate enters the constant 
TV in Eq.( 4.13| ) that is equal to the number of absorp- 
tion/emission events for the escape time and thus de- 
scribes the pumping strength. 

At T* — uj\[M E c the variance of persistent current 
fluctuations is stongly suppressed but still significantly 
depends on the bath temperature T even if T <C T* [see 
Fig.9]. 



<If> 




FIG. 9. Temperature dependence of the variance of the first 
harmonic of persistent current in units of (e/ro) 2 : the equi- 
librium temperature dependence (upper curve) and two close 
curves (blown up by a factor of 5) for the non-equilibrium 
case with T„ = 20E C corresponding to Q = 10, J\f = 4 and 
u> = 2, M = 100. The temperature dependence is significant 
in the non-equilibrium case even for T<T,. 

For instance at uj <C E c and T„ = uj^/M 3> E c we 
obtain 



dx 



-n 2 /{Ax) 



T 2 



sinh (irTx) 



(5.12) 



Eq.( |5.12 ) shows the dependence on the bath temperature 
T = TE C which is of the same type as in the absence 
of pumping Eq.(5.7), only the magnitude of fluctuations 
decreases by the factor of (E c /T*) 2 <C 1. So the overall 
width T* and the small steps in the electron energy dis- 
tribution of Fig. 5 manifest themselves in the variance of 
persistent current fluctuations. 



VI. DC CONDUCTANCE UNDER AC PUMPING 



size of the system L 3> Lt where Lt = y/D/T. The 
question we address here is whether or not the effective 
temperature T* is what should stand for the bath tem- 
perature T in the expression for Lt under ac pumping. 
The answer is not obvious in the geometry of an open 
quantum dot connected by the leads to electron reser- 
voirs [ see Fig. 10a, b]. The point is that there are several 
different electron energy distributions in such a problem. 
The electron energy distribution in each reservoir is sup- 
posed to be the equilibrium Fermi distribution with a 
certain chemical potential and temperature. In addition, 
there is the non-equilibrium electron energy distribution 
inside the dot under ac pumping. 

We define the Landauer conductance as the linear dc 
current response to the difference of the chemical poten- 
tials between two different reservoirs with the same tem- 
perature [Fig. 10a]. 

Another experimental situation corresponds to mea- 
surements of the linear dc current response to the per- 
turbation of the system 's Hamiltonian caused by the con- 
stant electric field inside the mesoscopic system. The 
corresponding response function will be referred to as 
the Kubo conductance. It can be realized as a current 
response in a ring that is pierced by a magnetic flux |3(J 
[see Fig. 10b]. 



f(— ) 

W 2T > 




(b) 



FIG. 10. Experimental geometries for the Landauer (a) and 
Kubo (b) conductances. In the case of Landauer conductance 
the dc current 7o is the response to the chemical potential 
difference pi — fi2 between two reservoirs which enter through 
the Fermi-distributions of the incoming electrons. In the case 
of Kubo conductance the dc current Io is the response to the 
perturbation jEot of the system's Hamiltonian. The passive 
lead connects the mesoscopic ring to the electron reservoir 
with the equilibrium distribution of electrons. 



In this section we consider the dc conductance in a 
mesoscopic system under ac pumping. This problem has 
been recently addressed in the work by Pedersen and 
Biittiker p8[ . Here we neglect the electron interaction 
and focus on the effect of heating by the ac field. It 
is well known that the mesoscopic conductance fluctua- 
tions are temperature dependent and decrease when the 



The flux is supposed to contain two parts: one is grow- 
ing linearly with time and causes the dc electric field. 
Another one produces the high-frequency pumping. 

We will show that these two cases are drastically dif- 
ferent (the same conclusion has been reached by Vavilov 
and Aleiner (29)). The Landauer conductance of non- 
interacting electrons is insensitive to heating and depends 
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only on the bath temperature T. At the same time the 
non-equilibrium electron energy distribution in the meso- 
scopic ring does matter for the Kubo conductance. In 
particular the mesoscopic fluctuations of the Kubo con- 
ductance should feel the effective temperature T» rather 
than the temperature of the electron reservoir T which 
the mesoscopic ring is connected to by the passive lead. 



A. Landauer conductance in terms of electron 
Green's functions in the time domain 

In order to prove the statement on the absence of sensi- 
tivity of the Landauer conductance to heating produced 
by ac pumping in the dot and to see the key difference 
between the Landauer and the Kubo conductances we 
re-derive the Landauer conductance to allow for an ar- 
bitrary ac pumping. To derive the expression for the 
Landauer conductance in terms of the electron Green's 
functions we proceed along the route used in Ref. O]. 



1. Formulation of the Landauer conductance 

For simplicity we consider a zero-dimensional dot de- 
scribed by the N x N random time- dependent matrix 
Hij(t) connected to two perfect semi- infinite (x < 0) 
leads each containing M channels labeled by a. There is 
neither disorder nor electron interaction with each other 
or with an external ac field inside the leads. Therefore 
the matrix Green's function G_ M j_ a a{t — t' , x — x') of the 



^ ± e ±ikx-iv F kt e ±i PF x^ xhen the wave function in the 
leads (x < 0) is: 



Keldysh technique Eq. ( 2£ ) for incoming electrons in the 
leads depends only on the difference of time and coor- 
dinates fjlf . Moreover, the incoming electrons are sup- 
posed to be in equilibrium at the bath temperature T and 
the chemical potential fi + 5/j,, where Sfi = ±{7/2 differs 
in sign for the leads 1 and 2. Then the Keldysh compo- 
nent G+ + a g(E, x — x') of the incoming electron Green's 
function G ++ a p(E, x — x') in the energy-coordinate rep- 
resentation takes the form [cf . Eq. (|2.9|)1 : 



G% +>a ,{E,x-x')=t^(^^ 
x [G* +>aP (E, x-x 1 )- G$ +>a(3 (E, x - x% 



(6.1) 



The retarded and advanced components for incoming 
electrons G^^Jt — t',x — x') in the leads are given 
by: 

G++, af} = i0(t - t') S aP 5(v F (t - t') - (x - x')) (6.2) 
Gi +M p = -i6(t' - t) 6 a p 6(v F (t - t') -{x- x')). (6.3) 

Here we linearized the Schrodinger equation near 
the Fermi momentum and introduced right (in- 
coming) and left (outgoing) movers *$>±\x,t) = 



* {a '>(x,t) = ^ ) {x,t) + ^>{x,t) 



(6.4) 



The Schrodinger equation for electron states inside the 
dot ipn{t) coupled to electron states in the leads ^^(x, t) 
is taken in the form: 

Y^liSnrndt - H nm (t)} 4> m (t) = * (Q) (0, t), (6.5) 



where W an is the 2M x TV coupling matrix. 

Neglecting electron-electron interaction we introduce 
the linear boundary conditions at x = 0: 



[* ( ?\0,t)-^(0,t)} = J2w am i> m (t), (6.6) 



For massive leads with the semiclassical electron motion, 
electrons adiabatically turn back at W am — * acquir- 
ing only a certain phase which may be included into the 
definition of ^ 



Equations Eqs.(6.5),(6.6) generate the corresponding 
equations for the matrix Green's functions: 

Y}i5 nm d t - H nm (t)]G ±:m/3 (t,t';x') (6.7) 

m 

= £ Wt a [G +± ^{t, f ; -x') + G_ ±>af3 (t, f ; -x% 



J2w am G ±tm/3 (t,t';x') = -iv F (6.8) 
m 

x [G + ± i0t /3(t, t'; —x) - G_± af3 (t, t'; —x)], 

where the label +(— ) corresponds to the incoming (out- 
going) electrons in the leads and G ± m/3 is the 'cross' 
Green's function of an incoming (outgoing) electron in a 
channel /3 in the leads and an electron at a site m in the 
dot. 

According to Eq.(^.Tj) the current I\,2{t) in the leads 1 
or 2 is given by the Keldysh components G+ + aa (t, t; 0) 
and G^_ aa (t, i;0) of the incoming and outgoing elec- 
trons: 

h,2 = iv F £ [G+ +aa {t,t;0)-G K ^ aa {t,t;0)}. (6.9) 



2. Analytical structure of the scattering matrix 



It is possible |ll|| to express the Keldysh component 
G K :_ a p(t, t'; 0) of the outgoing electrons in terms of the 
known Keldysh component G^ + a g(t, t'; 0) for the incom- 
ing electrons using the time-dependent scattering matrix 
S a p{t, t'): 
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G^ Q/3 (M';0)= (6.10) 
= J dh J dt 2 S ai {t,t l )G I l +n& {t 1 ^t 2 -Q)Sl fj {t 2 ,t'). 

It is crucial for us that the scattering matrix S a p(t,t') 
involves only the retarded component G^ n (t, t') of the 
electron Green's function inside the dot and the Keldysh 
component G^ m (t, t') drops out of the scattering matrix 

& 

SapW) = 8 aP S{t - - VF ttn G R m (t,t') Wl . 

(6.11) 

We remind that only the Keldysh component G^ m (t,t') 
contains an information on the electron energy distri- 
bution inside the dot. It describes the real transitions 
between energy levels in the dot which are constrained 
by the energy conservation law and the Pauli principle. 
In contrast to that the retarded component G R n (t,t') 
describes virtual transitions where no energy conserva- 
tion applies and thus the energy distribution function is 
irrelevant. 

Formally the fact that the scattering mat ri x is inde- 
pendent of G*„ (*,*') follows from Eqs.(|3),(U). Let 

us expr ess G |_ a g through Gj^_ n g an d G +mf3 using 

Eq.( |6.8|) and substitute it into Eq. j6.7j ). As the result 



of the transformation 2G 



appears in the r.h.s. of 
Eq. (|6.7|) and the matrix Hamiltonian for electrons in the 



++,a/3 



dot acquires an imaginary part H nm - 
H^ = H nm Tiv F 1 (W^W) r 



H n m , where: 



(6.12) 



According to the rules of the Keldysh technique the 
generic solution G+ mB to the transformed Eq. (p/fl) in- 
volves both the combination G^ n W^ a G^ + aj3 and the 
combination G^ ln W^ la G A + a3 , where G is the marix 
Green's function for electrons in the dot, for instance: 

G% m (t,t')=[id t -H R (t)]-l. (6.13) 

However, in this particular scattering problem we have: 

Gi +>a p{t,t';x-x' >Q)=0. (6.14) 

Eq.( |6.14 ) follows immediately from the S- function and 
the 0- function structure of Eq.(6.3). Therefore we ob- 
tain from Eqs.(|6.7|),(6.< 



G5+, a/ 9(M';0)= / dt 1 S a7 (t,t 1 )G£ fi7/3 (t 1 ,t , ;0). (6.15) 



with S a -y(t, t\) given by Eq.(6.11) which contains only the 
retarded component G^ TO of the electron Green's func- 
tion in the dot. 



In a similar way one can express G^_ a g(t, t'; 0) 
through G K : + a3 (t , t'; 0) using the Hermitean conjugated 
equations Eq.(|6.7|),(|6.8[) to finally arrive at Eq. fl6.10p , 



3. Expression for the Landauer conductance 



Now we substitute Eq.( |6.10D in Eq.fl6.9D with S a p(t, t') 
given by Eq.(6.11) to get for the linear dc response cur- 
rent I = h n (t): 

hn{t) = -V F 2 U { dh ( dt 2 [-»(<! - t 2 ) /(*! - t 2 )] (6.16) 



x Tr{(l + A)WG R (t,h)W^AWG A (t 2 ,t)W^} 

+ *(^V2) U J dh [-i(h - t) f(h - t)\ 

x Tr {(1 + A)W[G R (t,h) - G A (h,t)]W^} , 

where we have introduced the diagonal matrix A: 



-1, a € lead 1 
A„„ = i —1. a € lead 2 
0, otherwise 



(6.17) 



Equation Eq.(S.16) can be simplified using the identity: 

G R (t,t') - G A (t,t') = G A {[G-Y - [G-Y} G R 
= -2iVp X J dhG A (t,h)W^WG R (h,t') (6.18) 



which is obtained with the help of Eq.( 6.12 ). 

In order to obtain the dc Jo one has to average 
Eq.(EU6l) over t within the observation time T — > oo. 
This means an additional integration over t. Then using 
the fact that tf(t) is an even function one can r eplac e 
G A {t 1 ,t) by G A (t h) in the second term of Eq.fl6.16P. 
After that Eq.(6.18|) can be applied to yield: 



SLand — Vp 



dh / dn[-inf(n)) 



^ r l 2 dt_ 

-T/2 T 

x Tr {~W\\ + A) WG R (t, h)W*AWG A (h + t), t) 
+ W^WG R {t,h)W\l + A)WG A (h +77, t)} . (6.19) 

Finally we assume that the 2M x N coupling matrix W an 
has only 2M nonzero matrix elements, those with a — n; 



all these elements are taken equal to W. Then Eq.(3.1E) 
takes the form: 



ffLand — 



J dh J dr 1 [~ lV f\ V )}M- 2 Ti-{G R (t 7 h)G A (h+r 1 ,t) + G R (t,t 1 )G A (h+ V ,t)} 



(6.20) 
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where f(rj) is the Fourier transform of tanh(£'/2T); the 
overline denotes the average over time t; 7 = 2w^ 1 |W| 2 M 
is the electron escape rate, M is the number of channels 
in each lead and the subscripts 1 or 2 indicate that only 



sites connected to the first or the second lead should be 
taken into account in the summation over matrix indices. 



Using Eq.(6.18) one can |4j recast Eq.(6.20) as a sum 
of two parts ^Land = <?i + 92 where 



.91 = ^7 / dr) [-ir) f(r))] M- 1 tr {c«(t + r,, t) - G A (t + r,,t)} , 



and 



.92 



\l 2 J dt x J dr, [-i v f(rj)] M- 2 tr { AG*(t, t 1 )AG A (t 1 + v ,t)} 



(6.21) 



(6.22) 



In Eqs.( 6.21 ),( 6.22 ) the symbol tr denotes the matrix 
trace taken onl y ove r sites connecte d to leads and the 
matrix A in Eq.( 6.22 ) is given by Eq.( 6.17 ). 

The first contribution g\ is proportional to the local 
density of states in the regions of leads. It describes the 
effect of electron escape in each lead separately. The sec- 
ond term g 2 describes the mutual effect of both leads. It 
is similar to the conventional term jG R jG A in the Kubo 
conductance, since Trj = trA = 0. 



B. Why there are no loose diffusons in the problem 
of Landauer conductance under ac pumping 

We stress once again that the Landauer conductance 
is independent of the Keldysh component of the electron 
Green's function in the dot. This is already an indication 
that elecron kinetics inside the dot under ac pumping is 
irrelevant for the Landauer conductance. In particular 
this means that the electron energy distribution function 



f{rf) i n Eq. (6.20) is not renormalized in the way similar 
to Eq.(fh~L2l). 



Formally this follows from the difference in the struc- 
ture of the retarded- advanced junctions. Let us com- 



pare the retarded-advanced junction Eq.(3.5) and that 
in Eq.( 6.22 ). The difference is that Eq^E^fj) contains the 
electron interaction with the ac pumping field rdt 1 A(ti) — 
r<9 tl A(ii + rj) or jA(ti) — jA(ti + ry). At the same 
time Eq.( 6.22 ) contains — ir]A oc A[A (£i) — Ao(ti + 
rj)] which stems from the interaction with the dc field 
Ao(t) = — Eo< that causes the chemical potential differ- 
ence U. Correspondingly, the expression for the loose 
diffuson analogous to Eq.(4.6) is proportional j32j to 



A p o = Em^W^m^i" 1 ) oc trA = (instead of 
r^o)- An expansion of G R,A in powers of the ac pump- 
ing field analogous to that leading to the triangle in Fig. 4 
does not help. Because the retarded-advanced junction 
in Eq.( |6.22 ) does not contain the ac field the correspond- 
ing triangle is linear in the ac field and vanishes after 



averaging over time t [cf. E q.(4.11 )]. Thus the loose 
diffuson corresponding to Eq.( 3.22 ) vanishes in the same 
way as the one that corresponds to the conventional term 
jG R jG A in the Kubo conductance. The structure of 
Eq.( |6.2l|) does not allow to build the loose diffuson be- 
cause it does not contain the retarded-advanced junction 
al all. 

Since the heating effect is associated with the loose 
diffusons the above arguments allow one to conclude |2£( | 
that the Landauer conductance of non-interacting elec- 
trons is not sensitive to heating. 



C. Sensitivity of the Kubo conductance to heating 

However, this is not true for the Kubo conductance 
which is sensitive to heating. 

Consider an experimental situation shown in Fig. 10b. 
Here the electric field Eq is produced by the linear in time 
component of the magnetic flux piercing the mesoscopic 
ring 4>o — LAo = —LEgt. The ac pumping is produced by 
the oscillating part of the magnetic flux <p ac {t) = A(t)L. 
In the framework of the dynamical approach the steady- 
state regime is only reachable if the mesoscopic system 
is connected to the reservoir by a passive lead. It sup- 
ports no net current but results in the particle exchange 
between the mesoscopic system and the reservoir neces- 
sary to remove heat produced by the ac pumping. Phe- 
nomenologically, the effect of the passive lead is described 
by the escape rate 7. This geometry corresponds exactly 
to the formalism developed in Sec. II and Sec. Ill, where 



H c - { (t) = Ha + Hi =jEoi-jA(i). 



(6.23) 



Now each ray in Fig. 2 represents a sum Tlo + Ti\. Since 
we are interested in the linear response to Tlo, one ray 
is special: it corresponds to Tlo while all other rays cor- 
respond to Hi- Then using Eqs.(2.8)-(3.6) one obtains 
Iq = hn(t) = [at + a 2 + 0-3] E , where: 



0-1 



d»?/(»7)Tl{jGJ 



■(t, ii)j [-tti]G«(ti, t - 77) - jG A (t + 77, ti)j [-it! 



G A (h,t)} 



(6.24) 
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(T2 = dt 



dr, [-ivKv)} Tr {jG*(i, hfiG^h + V ,t)} , 



<7 3 = / ^ 1 Tr{jGW(t,t 1 )jHii]G A (ti^)+jG fl (i,<i)jHii]G (a) (^i^)}- 
Here G fi, ^(i,i') are Green's functions in the mesoscopic system subject to ac pumping and 
G<°»(M')= / dt" / ^/>y)G fl (t,i" + r;)[H 1 (t" + r / )-H 1 (i")]G' 4 (t",t'). 



(6.25) 
(6.26) 

(6.27) 



The first two parts Eqs. fl6.24 ), ( 6.25 ) do not allow to 
build the loose diffuson for exactly the same reasons as 
Eqs.(3.21),( |6.22 ) for the Landauer conductance. How- 
ever, the 'anomalous' Green 's function in Eq.(|6~2^) 
contains the retarded-advanced junction of the same 
structure as we discussed in Sec. IV which does allow to 
build the loose diffuson. This is the part where the heat- 
ing effects oridinate from. 

One can see the direct analogy with the situation dis- 
cussed in Sec.Vb. The role played by the current I^ 2 \ 
Eq.( |5.9| ), is now played by the part of conductance 173. 
The parent diagrams for the problem of conductance fluc- 
tuations are the usual two-diffuson and two-cooperon dia- 
grams considered in Refs. p^ j. All the daughter-diagrams 
with one or two additional loose diffusons can be ob- 
tained from the disorder averages {(Ji^gs) and (of), re- 
spectively. It is not difficult to show that all of them 
again lead to the renormalization of the energy distribu- 
tion function given by Eq.(4.12). 



D. The difference between the Landauer and the 
Kubo conductance 

Let us discuss the difference between the experimental 
situations described by the Landauer and the Kubo con- 
ductances. The principal difference between them is that 
the quantity we call the Landauer conductance is the lin- 
ear current response to the variations of the parameters 
of the electron reservoirs (the chemical potential differ- 
ence) while the Kubo conductance is the linear current 
response to the variation of the system's Hamiltonian 
(the term 7io in Eq.( 3.23 )). We stress that the difference 
lies in the physical situations and not in the methods of 
description. For instance, the photovoltaic effect can be 
considered both in the framework of the scattering ma- 
trix approach JHJ and in the nonlinear response approach 
similar to the one we used in Sec.V for the persistent cur- 
rent fluctuations. In both cases the result is expressed 
in terms of the renormalized electron energy distribution 
Eq.(4.12). The reason is that the photovoltaic current 
is a non-linear response to the variation of the system's 



Hamitlonian (fH\ in Eq.(3.23)) caused by the ac pumping 
field which should be treated in the same way as the vari- 



ation TLq in Eq. ( |6.23|) that causes the Kubo conductance. 
That is why in terms of the sensitivity to heating the pho- 
tovoltaic effect is similar to the Kubo conductance and 
drastically different from the Landauer conductance. 

We note that the key difference between these two 
situations examplified by Fig. 10a and Fig. 10b is that 
the initial density matrix at t = — 00 prior to switch- 
ing on the interaction with external ac field (which is a 
starting point in the Keldysh formalism |l5}]), is already 
non- equilibrium in the problem of Landauer conductance 
[Fig. 10a], since it involves two Fermi-distributions with 
different chemical potentials. At the same time in the 
situation represented by Fig. 10b the initial density ma- 
trix is equilibrium. 



VII. CONCLUSION 

In the present paper we have considered a general for- 
malism of the dynamical approach to nonlinear response 
of mesoscopic systems. For completeness of presentation 
and for tutorial purposes we have shown how causality 
principle leads to the powerful machinery of analytical 
continuation |l2| which is most conveniently realized in 
the matrix algebra of the Keldysh technique P^ , [15| , |liT| . 
The principal goal of the paper was to demonstrate the 
capability of the dynamical approach and its apparent 
limitations. We have shown that the diagrammatic tech- 
nique developed in Ref. f2^,|2j| automatically describes 
the electron diffusion in the energy space under the ac- 
tion of ac pumping. We considered few simple exam- 
ples of how the renormalization of the electron energy 
distribution occurs due to the loose diffusons and have 
demonstrated its equivalence to the solution to the ki- 
netic equation. 

We have shown that the quadratic response in the 
closed mesoscopic systems is singular in the dynamical 
approach because of the singularity of the loose diffusons 
in the quadratic in the pump field approximation. Yet 
in contrast to Ref. |l8| this does not lead to an infinite 
dc current arising under ac pumping, as the quadratic in 
A approximation breaks down well before the singular- 
ity develops itself. The situation is somewhat opposite 
to that described in Ref. |Tq]. The ensemble-averaged 
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dc current does not depend on the singular loose dif- 
fusons whatsoever Jl9|] and thus is finite anyway. The 
mesoscopic fluctuations of the dc current in an isolated 
ring are zero because the corresponding effective temper- 
ature T* — > oo in the absence of dissipation or an elec- 
tron escape. This clearly shows that the dynamical ap- 
proach based on the 'minimal model' of non-interacting 
electrons in an impurity potential intercating with the 
external classical field is insufficient for describing the 
closed systems. 

However, this model is reasonable for open mesoscopic 
systems connected with the electron reservoir by massive 
leads. We have carefully studied the case of the Landauer 
conductance in a quantum dot under ac pumping, since 
here one can see the danger of ad hoc replacement of the 
bare electron energy distribution by the renormalized one 
with the simultaneous deletion of all the diagrams with 
the loose diffusons. The correct diagrammatic analysis 
shows that the diagrams with the loose diffusons do not 
arise in this problem, and the electron energy distribution 
that enter e.g. the variance of conductance fluctuations 
stays unrenormalized by heating. At the same time, the 
same type of diagrammatic analysis shows that the vari- 
ance of the Kubo conductance is sensitive to the renor- 
malization of electron energy distribution. This sets yet 
another borderline between these two formulations of the 
problem of conductance and shows their deep physical 
difference and significance of the experimental geometry 
for conductance measurements. 

Though the results presented in the paper are purely 
perturbative, there is a bridge to non-perturbative 
schemes. The obvious extension of the present theory 
is the time-dependent random matrix theory (TRMT) 
which applies to the problem of quantum dot under ac 
pumping. We found that in contrast to the equilib- 
rium or the linear response th eories , there a re tw o dif- 
ferent TRMT formulations Eq. fl3.37D and Eq. fl3.32|) , one 
for the adiabatic pumping with the typical frequency uj 
smaller than the Thouless energy E c and another one for 
the case of high frequencies to 3> E c . The correspond- 
ing matrix Hamiltonians have different symmetry: they 
are real symmetric for the adiabatic case and contain 
an imaginary anti-symmetric time-dependent part (with 
zero time-average) in the case of high frequencies. 

Another obvious extension is the functional formula- 
tion in terms of the non-linear sigma-model In this 
connection we note that in all field theories with the cor- 
rect vacuum, the loose propagators may not arise. In 
the particular formulation of the non- linear sigma-model 
based on the Keldysh formalism the electron energy dis- 
tribution is a part of the vacuum solution. However, the 
example of the Landauer conductance (that depends on 
the bare energy distribution function ) and the Kubo 
conductance or the photovoltaic effect (that depend 
on the renormalized energy distribution function) shows 
that the vacuum could be non-unique or containing dif- 



ferent sectors which enter in a different way in one or 
another observable. 

Finally we mention the effect of electron-electron in- 
teraction. It seems plausible that in the presence of an 
electron-electron interaction the Landauer conductance 
will be sensitive to heating. Another effect of electron 
interaction is the inelastic electron scattering. In our dy- 
namical approach we have neglected this effect at all. 
In open mesoscopic systems this has resulted in the elec- 
tron energy distribution Eq.(4.12) of a very peculiar form. 
The characteristic feature of this distribution is that it 
depends not only on the effective temperature T* but also 
on the bath temperature T « T,. The consequence of 
this two-parameter dependence takes its extreme form 
in the temperature dependence of the persistent cur- 
rent or phot ovoltaic current fluctuations. In particu- 
lar, Eq.( |5.12 ) shows how the separation of parameters 
occurs in the variance of the persistent current fluctua- 
tions under ac pumping. The parameter T* turns out 
to determine only the pre-factor in front of the function 
that depends on the bath temperature T. Such a be- 
havior is only possible if electron-electron interaction is 
neglected. Inelastic processes due to electron-electron in- 
teraction are always favorable to the Fermi-distribution 
function with some effective temperature T e i . They work 
to diminish all the deviations from the Fermi-distribution 
and should certainly reduce the dependence of the persis- 
tent current fluctuations on the bath temperature. This 
(or similar temperature dependences for the photovoltaic 
effect) can be a useful tool for the experimental investiga- 
tion of electron-electron interaction in an open quantum 
dot. 
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